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Abstract 


The combination of trellis coding and multiple phase-shift-keyed (MPSK) 
signaling with the addition of asymmetry to the signal set is discussed with 
regard to its suitability as a modulation/coding scheme for the fading mobile 
satellite channel. For MPSK, introducing nonuniformity (asymmetry) into the 
spacing between signal points in the constellation buys a further improvement 
in performance over that achievable with trellis coded symmetric MPSK, all 
this without increasing the average or peak power, or changing the bandwidth 
constraints imposed on the system. 

Whereas previous contributions have considered the performance of trellis 
coded modulation transmitted over an additive white Gaussian noise (AWGN) 
channel, the emphasis in this paper is on the performance of trellis coded 
MPSK in the fading environment. The results will be obtained by using a 
combination of analysis and simulation. It will be assumed that the effect of 
the fading on the phase of the received signal is fully compensated for either 
by tracking it with some form of phase-locked loop or with pilot tone 
calibration techniques. Thus, our results will reflect only the degradation 
due to the effect of the fading on the amplitude of the received signal. 

Also, we shall consider only the case where interleaving/deinterleaving is 
employed to further combat the fading. This allows for considerable 
simplification of the analysis and is of great practical interest. Finally, 
the impact of the availability of channel state information on average bit 
error probability performance is assessed. 
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SECTION I 
INTRODUCTION 


There is a growing need for reliable transmission of high-quality voice 
and digital data in satellite-based land mobile communication systems. These 
systems, which will be part of an emerging all-digital network, are both power 
and bandwidth limited. To satisfy the bandwidth limitation, one can employ 
bandwidth-efficient modulation techniques such as those that have been 
developed over the past several years for terrestrial microwave communications 
systems. Examples of these are multiple phase-shift keying (MPSK), quadrature 
amplitude modulation (QAM) and the various forms of continuous phase frequency 
modulation (CPM). When power is limited, forward error correction (FEC) 
coding is ordinarily used. 

When both limitations are imposed simultaneously, as in the mobile 
satellite application, it is most often not possible to achieve the desired 
throughput with either technique acting alone. Instead, what is required is 
the integration of a bandwidth-efficient modulation scheme with some form of 
FEC coding to exploit the best possible attributes of both. 

In the past, coding and modulation were treated as separate operations 
with regard to overall system design. In particular, most earlier works on 
coded digital communication systems independently optimized: (a) conventional 

(block or convolutional) coding with maximized minimum Hamming distance, and 
(b) conventional modulation with maximally separated signals. 

About a decade ago, using random coding bound arguments, it was shown 
that considerable performance improvement could be obtained by treating coding 
and modulation as a single entity [1]. Many years later, this concept was 
formalized into a rigorous theory which showed that optimally designed rate 
n/(n+l) trellis codes suitably mapped (to maximize Euclidean distance) into 
conventional 2 n+ ^-point signal sets can provide significant coding gain 
without bandwidth expansion when compared with uncoded conventional 2 n -point 
signal sets [2]. It is this work that has laid the foundation for the design 
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and development of all power and bandwidth-efficient digital modems found in 
practice today and those that are to come in the future. 

The most common application of such trellis coded modulation (TCM) 
techniques is in the new generation of modems being developed for the 
telephone channel. Indeed, the present state of the art. is a rate 6/7, 

8-state trellis coded 128-point QAM which is capable of transmitting 14.4 kbps 
(2.4 ksps) over good-quality (Dl-conditioned or better) leased telephone 
lines [3]. 

Thus, if it is practical to send 14.4 kbps over the telephone channel, 
transmitting 4.8 or 9.6 kbps information over a 5-kHz satellite channel 
(typical of present considerations) might appear to be simple. 

Several reasons make this supposition untrue, most of which relate to the 
additional sources of degradation present on the mobile satellite channel but 
absent on the telephone channel. First, Doppler frequency shifts due to 
mobile vehicle motion can be a serious source of performance degradation if 
not compensated for. Second, the fact that the 5-kHz mobile channel is 
actually a slot in a frequency-division multiple access (FDMA) system brings 
on the problem of interference due to energy spillover from adjacent 
channels. Third, the satellite channel is inherently a nonlinear one 
primarily due to the HPA in the transmitter. Thus, one must either employ 
constant envelope modulations and operate at full power or, if using a 
nonconstant envelope, but bandwidth-efficient modulations such as QAM, the HPA 
operating point must be backed off in power to produce an approximately linear 
channel. The most serious source of impairment that does not exist on the 
telephone channel is the combination of multipath fading and shadowing, i.e., 
for reliable performance, the system must combat short fades and recover 
quickly from long fades. Fading, which for mobile satellite channels can be 
assumed to be modelled by a Rician distribution, not only introduces an error 
floor into the system but also makes the problem of carrier recovery more 
difficult. Depending on the ratio of direct and specular (coherent component) 
to diffuse (noncoherent component) signal power, one might even be required to 
employ differentially coherent or noncoherent detection techniques, thus 
sacrificing the power saving associated with coherent detection. Finally, 
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even if the above sources of degradation were absent, the power limitation 
imposed by the mobile satellite channel would preclude transmission at the 
high data rates achievable on the telephone channel. 

Also, whatever technique is decided upon must be able to be implemented 
and installed in the vehicle with a minimum of cost and complexity, perhaps 
two orders of magnitude less than that associated with a telephone channel 
modem. 

This report shows that the combination of trellis coding and MPSK 
signaling with the addition of asymmetry to the signal set is a desirable 
modulation/coding scheme for the fading mobile satellite channel. For MPSK, 
introducing nonuniformity (asymmetry) into the spacing between signal points 
in the constellation (see Figure 1 for example) buys a further improvement in 
performance over that achievable with trellis coded symmetric MPSK, all this 
without increasing the average or peak power, or changing the bandwidth 
constraints imposed on the system [4,5]. 
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Figure 1. Symmetric and Asymmetric MPSK Signal Sets 
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Whereas previous contributions [2-5] have considered the performance of 
trellis coded modulation transmitted over an additive white Gaussian noise 
(AWGN) channel, the emphasis in this report is on the performance of trellis 
coded MPSK in the fading environment. The results will be obtained by using a 
combination of analysis and simulation. It will be assumed that the effect of 
the fading on the phase of the received signal is fully compensated for either 
by tracking it with some form of phase-locked loop or with pilot tone 
calibration techniques [6,7]. Thus, our results will reflect only the 
degradation due to the effect of the fading on the amplitude of the received 
signal. Also, we shall consider only the case where interleaving/ 
deinterleaving is employed to further combat the fading. This allows for 
considerable simplification of the analysis and is of great practical interest. 
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SECTION II 
SYSTEM MODEL 


Figure 2 is a block diagram of the end-to-end system under 
investigation. Input bits representing data or digitally encoded speech are 
passed through a rate n/(n+l) trellis encoder which is normally implemented 
with a combination of n shift registers (the memory of the encoding operation) 
and appropriate mod-2 adders (XOR gates). The encoder output symbols are then 
block interleaved to break up burst errors caused by amplitude fades of 
duration greater than one symbol time. While in practice the depth of 
interleaving is finite and chosen in relation to the maximum fade duration 
anticipated, for the purpose of analysis, we shall make the usual assumption 
of infinite interleaving depth. This assumption provides a memoryless channel 
for which well-known bit error probability bounding techniques can be 
applied. Our simulation results will, however, reflect a finite interleaving 
depth. Thus, these results will be slightly pessimistic when compared with 
those derived from theory. The primary purpose of the analysis is to indicate 
through simple example the trend of the performance behavior as a function of 
the various system parameters leaving the actual numerical performance to be 
predicted by the software simulations. 

Groups of n+1 interleaved symbols are mapped (with a ROM) into the MPSK 
signal set according to the set partitioning method [2]. The in-phase and 
quadrature components of the mapped signal point are digitally pulse shaped 
[to limit intersymbol interference ( IS I ) ] and modulated onto quadrature 
carriers for transmission over the channel. If pilot tone calibration 
techniques [6,7] are used to recover the faded carrier at the receiver, then 
the pilot tone (or tones) must be added to the above data-modulated signal 
before transmission. 

At the receiver, the faded, noise-corrupted in-phase and quadrature 
signal components are demodulated with the extracted pilot tone(s), q-bit 
quantized for soft decision decoding, and then block deinterleaved . The 
metric chosen for the Viterbi algorithm in the decoder depends upon whether or 
not channel state information (CSI) is provided [8]. As indicated in 
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PILOT TONE(S) 
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Figure 2. System Block Diagram 


















Figure 2, a measure of CSI can be obtained from the power in the recovered 
pilot tone(s). Furthermore, the number of bits of quantization, p, for this 
operation can be much smaller than q since the accuracy of the CSI has only a 
secondary effect when compared with that of the soft decisions themselves. 
Finally, the tentative soft decisions from the Viterbi decoder are stored in a 
buffer whose size is, in our application, a design parameter. In particular, 
for the case of speech transmission, the total coding/decoding delay must be 
kept below about 60 ms so as not to be objectionable to the listener. Thus, 
for a given input bit rate, the decoder buffer and interleaving frame sizes 
must be limited so as to' produce at most a 60-ms delay. Again for simplicity, 
we shall assume an infinite buffer in the analysis, whereas the simulations 
will reflect a finite buffer size in accordance with the delay limitation. 
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SECTION III 

DERIVATION OF THE PAIRWISE ERROR PROBABILITY BOUND 


The basic analysis model for the system of Figure 2 is illustrated in 
Figure 3. The box labelled encoder is actually the combination of the trellis 
encoder (with binary input and output) and the mapping function of .Figure 2 
which, as shown in [9] can be described by a trellis encoder with binary input 
and M-ary output. In our notation, M = 2 n+ ^. Also, if the input bit rate 
is denoted by R^, then the symbol rate input to the channel is R^/n. We 
denote a coded symbol sequence of length N by 



( 1 ) 


where the kth element of x, namely x^, represents the transmitted MPSK 
symbol at time k and is a nonlinear function of the state of the encoder s^ 
and the n information bits, u, , at its input, vis.. 


\ - f( V Slc> 


( 2 ) 


The next state of the encoder 
the present state s fc and the n 
described by 


s, , at time k+1 is a nonlinear function of 
k+1 

-bit input u, which is mathematically 


s k + i * 8( v V 


(3) 


Corresponding to x, the channel outputs the sequence 

X - < y i. y 2 ’ V 


(4) 



Figure 3. General Memoryless Channel 
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where the kth element y^, representing the output at time k, is given by 

y k - p k*k * % <5) 

In (5), is a normalized (unit mean-squared value) random amplitude 

whose probability statistics depend on the fading channel model (to be 

discussed later on) and n, is a sample of a zero-mean Gaussian noise process 
2 * 

with variance a . If side information is available, the corresponding 
side information sequence z will be denoted by 

z = (z., z 2 , . . . , z N ) (6) 

As stated in the introduction, we shall assume adequate (theoretically 
infinite) interleaving and deinterleaving so that the coding channel is 
memoryless. Under this assumption, the p^'s are independent random 
variables (r.v.'s).and hence the channel probabilities satisfy 


P N (*|x.z) -n^n'V’n’ 
n=l 


and 


N 


q„<5> =|~[ Q(« n > 

n=l 


(7) 


( 8 ) 


For any coded communication system, the decoding process uses a metric of 
the form m(y,x; z) if side information is available and m(y,x) if it is not. 
Whatever metric is selected, it is desirable from the standpoint of 
simplifying the decoding process that it have an additive property, namely 
that the total metric for a sequence of symbols is the sum of the metrics for 
each channel input and output pair. In terms of (1), (4), and (6), this takes 
the form 


N 


m(y,x; z) = ^ ] 


m(y 


Z n ) 


n=l 


(9) 
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The maximum-likelihood metric 
m(y,x; z) = fi-n P N (y|x,z) 


( 10 a) 


when side information is available or 
m(y,x) = £n p N (y|x) 


... (10b) 


when no side information is available satisfies the requirement in (9). This 
is easily seen by substituting (7) in (10) and recalling that the logarithm of 
a product equals the sum of the logarithms- Thus, we shall use this metric in 
what follows. Also for simplicity, we shall use only the notation of (10a) 
keeping in mind that the case where no side information is available is a 
special case where the metric m(y,x; z) does not depend on z. 

To find an upper bound on the average bit error probability performance of 
the Viterbi decoder, we must first find the pairwise error probability 
P(x ■* x) which represents the probability of choosing the coded sequence 

A A A A 

x = (x]_, X£, Xfl) instead of x = (x]_, X 2 , ...» xpj). Since the decoder 

incorrectly decides x is the transmitted sequence when 

m(y,x; z) > m(y,x; z) ( 11 a) 

or, equivalently 


N 


N 


z ) m(y ,x ; z ), 

n n n n n n 


n=l 


n=l 


(lib) 


the pairwise probability P(x -* x) is given by 


N 


N 


P(x x) = Pr < 7 m(y ,x ; z ) > \ 1 m(y ,x ; z ) lx 1 
_ — \ n n . n - n n n ' — j 


n=l 


n=l 


( 12 ) 
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Applying the Chernoff bound, (12) becomes [10] 

E{exp (X[m(y n ,x n ; z fl ) - m(y n> x n ; z n >])|x} (13) 

nen 

where "E" is the statistical expectation operation, X is the Chernoff 
parameter to be optimized, and n is the set of all n such that ^ x^. 

To simplify (13) further, we must specify whether or not side information is 
available. 


P(x * x) < 


n 


A. IDEAL CHANNEL STATE INFORMATION 


The assumption of ideal channel state information is tantamount to assum- 
ing that the side information r.v. z^ is equal to the fading amplitude p^. 

Thus, since n is Gaussian distributed and in general, x and y are 
n e n J n 

complex quantities, then using (5) and (10a) gives 




(U) 


2 

where for simplicity we shall ignore the 1/2 a factor since it would later be 
absorbed in the Chernoff bound parameter X anyway. Substituting (14) into (13) 
and conditioning on [p we get 


P(x ■* x|g) < | [ exp [-Xp 2 |x - x | 2 ] E{exp [-2Xp Re{n (x - x )*}]} (15) 

— — — l l r r n n n n-nn n 

nen 

Representing the complex noise r.v. n^ in terms of its real and imaginary 

parts where Re(n } and Im{n } are uncorrelated zero-mean Gaussian r.v.'s, 
n 2 2 n 

each with variance o„ = a / 2, it can be shown that 

N 


E{exp [-2Xp Refn (x - x )*}]} = exp [2X 2 p 2 a 2 |x - x | 2 ] 

1 r n 1 n n n JJ r r nN'n n 1 

Substituting (16) into (15) gives after some simplification 

P(x ♦ x|^) < J~T exp (-Xp^!x n - x n | 2 (l - 2Xo 2 )] 


nen 


(16) 


(17) 
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2 •' * 

Since 0 ^ is independent of n, optimizing (17) over the Chernoff bound 

parameter X yields 
1 

\>pt “ 4c 2 (18) 

which when substituted in (17) produces the desired result, namely, 

p 2 |x - x | 2 

P(x -*• x|^) < | | D n n n ; D = exp (-l/8o 2 ) (19) 

nep 

Finally, (19) can be rewritten as 

P(x -» x|^) < D d2 (?»?) ( 20 ) 

where 

a 2 (x,i)^ pX-ij* < 21 > 

nen 

represents the square of the weighted Euclidean distance between the two 
symbol sequences x and x. 

Finally, the pairwise error probability upperbound is obtained by 
averaging (20) over the probability density function (p.d.f.) of 

B. NO CHANNEL STATE INFORMATION 

When no channel state information is available, then the metric of (10b) 
becomes 


m ( y n ,K n ) = - 




( 22 ) 
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Substituting (5) into (22), then, analogous to (15), we now get 


P(x -► x|£) < J~ J exp [-X|x n - x n | 2 - 2X(p n - 1) Re{x n (x n - x n )*}] 


nen 


x Efexp [-2Xp Re{n (x - x )*}]} 
1 r n 1 n n n J 


( 23 ) 


Again using (16), Eq. (23) simplifies to 


P (x + x 


le.) < Y\ exp t_X l X n “ X J 2(1 " 2Xa N ) 


nen 


- 2X(p - 1) Re{x (x - x )*}] 
n n n n 


(2M 


2 -> 2 

For constant envelope signal sets such as MPSK where |x| = |x| , (2U) can be 

further simplified by noting that 


| x - x| 2 = 2Re{x(x - x)*} 


(25) 


Thus, substituting (25) into (24) and renormalizing the Chernoff parameter 

2 

(i.e., replace X by 2X0^), we get the desired result analogous to (17), 
namely. 


Ms ■> slft,X> < D c 2 < S'Slfi-' X) (26) 

where D is again given by (19) and 

c 2 (x,x|£,X) = V' 4X(p - X) | x - x | 2 
— ~ n n n 

p lx - x I 2 - 4X 2 d 2 (x,x) (27) 

(26) cannot be optimized over X to yield a constant 
Thus, in this case, we must first average over the 

fading distribution. 


nen 

nen 

Note that unlike (17), Eq. 
value for this parameter. 
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SECTION IV 

DERIVATION OF THE BIT ERROR PROBABILITY BOUND 

To derive the upper bound on bit error probability from the pairwise 
error probability bound, we follow the transfer function approach taken in 
[5]. In particular, we first find the unconditioned pairwise error 
probability by averaging (20) or (26) over the p.d.f. of p. (This averaging 
will be denoted by an overbar.) When this is done, the pairwise error 
probability can be expressed in the form* 


P(x -> x) = D A 

where for the ideal channel state information case 

A = \ |x - x | ^ 

/ j r n' n n 1 

nen 

and for the no channel state information case 



nen 


(28) 


(29a) 


(29b) 


In dealing with upper bounds of this type, it is convenient to work with 
a pair-state transition diagram [11] where the pair-state and pair- 
information symbol are defined as 


S 



\ £ <v V 


(30) 


*For uniformity of notation, we shall drop the dependence of the pairwise 
error probabilty on X for the no channel state information case. 
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where and are the estimates for the state of the decoder and the 
information symbol, respectively. We are in a correct pair-state when s^ = s^ 
and in an incorrect pair-state when s^ ^ s^. Substituting (2) into (29), we 
can write the latter as 

A = V' 6 2 (S ,U ) (31) 

f -j n n 

nen 


where 

6 2 (S ,U ) = p 2 |f(s , u ) - f(s , u )| 2 (32a) 

nn n nn nn 

or 

S 2 (S ,U ) = 4X(p - X) | f (s ,u ) - f ( s ,u )| 2 (32b) 

n n r n n n n n 

When no fading is present, (28) is given by [5] 

P(x -► x) = D a (33) 

where A is still of the form in (31) with, however, 

6 2 (S ,U ) = | f (s ,u ) - f(s ,u )| 2 (34) 

n n 1 n n n n 1 

In terms of the above definitions, it can be shown by analogy with the 
results in [5] that the average bit error probability is upper bounded by 


P b 1 n di ?<D,I) ll=l 


(35) 


where T(D,I) is the transfer function of the pair-state transition diagram 
whose branch label gains are modified from those for the no fading case as 
follows. In the absence of fading each branch label gain has a factor 
(S n ,U n ) where S 2 (S n ,U n ) is given by (34). For the fading case, we 
simply replace D^ 2 (S n ,U n ) by D^ 2 (S n ,U n ) where again the overbar denotes 
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2 

averaging over the p.d.f. of p and 6 (S n ,U n ) * s 8^ ven ^ (32a) or (32b) as 
appropriate depending on the availability of channel state information. 

Finally, for the case of no channel state information, we must minimize 
the upper bound of (35) over the Chernoff parameter to obtain the tightest 
upper bound. Recall that for the ideal channel state information case, we 
were able to perform this optimization at the pairwise error probability level. 
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SECTION V 

CHARACTERIZATION OF THE FADING CHANNEL 


For mobile satellite communication, multipath fading produces a received 
signal with an amplitude* which can be modelled by Rician statistics with 
parameter K representing the ratio of the power in the direct (line-of-sight) 
and specular component to that in the diffuse component. If shadowing is 
severe, then a Rayleigh statistical model becomes appropriate which can be 
looked upon as the limiting case of a Rician channel when K approaches zero. 
Of course, the case of no fading corresponds to a Rician channel with K 
approaching infinity. 

Mathematically speaking, the above statements correspond to a p.d.f. for 
the fading r.v. given by 


( 2p(l + K) exp [-K - p 2 (l + K)]I n (2p yK(l + K)); p > 0 
p(p) = (36) 

(0 ; otherwise 


where Iq(x) is the zero-order modified Bessel function of the first kind. 


*We assume that the phase distortion produced by the fading is fully compen- 
sated for either by tracking it with some form of phase-locked loop or with 
pilot tone calibration techniques [6,7]. 
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SECTION VI 
AN EXAMPLE 


Consider the case of rate 1/2 trellis coded asymmetric QPSK using a 
2-state trellis. The appropriate set partitioning is illustrated in Figure 4, 
the trellis diagram in Figure 5, and the pair-state transition diagram in 
Figure 6. The performance of this system in the absence of fading was treated 
in [5] with the following results: 


T(D, I ) 


4ac 

1 - 2b 





(37a) 


or 


T(D, I ) 


ID 


4(l+2a) 

1-KX 



(37b) 


where D as defined in (19) becomes [5] 


D 



— = system bit energy-to-noise spectral 
0 density ratio 


(38) 


and a is the ratio of powers between the I and Q channels which is related 
to the angle <j> that defines the asymmetry (see Figure 4) by 


a = 



(39) 


Substituting (37) into (35) gives 
4(l+2a) 



( 40 ) 
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Optimizing (AO) over the asymmetry produces 


4 to D _ V N 0 
a ■ an 3 an 3 

which when substituted into (40) gives the desired upper bound 
P b < ^ exp (-2^/^) 


( 41 ) 


(42) 


A. IDEAL CHANNEL STATE INFORMATION 


Recalling (29a), we see that the transfer function T(D,I) of the 
pair-state transition diagram for the case of ideal channel state information 
is obtained merely by replacing by d p 2 b 

in the branch label gains of (37a) 

where B = 4, 4/(1 + a), or 4a/(l + a) as appropriate. For the Rician p.d.f. of 
(36), DP 2 * 3 evaluates to* 


— r / v -g- K -.. r 

D p fi _ / I tl .. n 1 +K+B Y . _b_ 

° " \1 + K + B Y y u ’ Y " 4 N q 


which for the Rayleigh special case (K = 0) becomes 


(43a) 


D 


2 

pi _ 1 

1 + By 


(43b) 


Evaluating T(D,I) of (37a) using (43) and performing the differentiation 
required in (35) gives the upper bound on as 


P b i 


¥ 3 ° 


*i * 3 



( 44 ) 


*Note that E^ represents the total average received power (direct plus 
specular plus diffuse) in the data bandwidth. In the general case of an 
n/(n+l) code rate, E^ in the definition of X would be replaced by nE^. 
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where 


r *i K . - 1 +K 

M 1 + K + G.y ’ M ~ 1 + K + B.y 


; i s 1* 2, 3 


lx 

Q, = 4 ; S. = t-=— ; fl. = ~— 
1 2 1 + a 3 1 + a 


( 45 ) 


To obtain the best performance in the presence of fading, one should 
optimize (44) over the asymmetry parameter a. Before doing this, however, we 
shall first examine the behavior of (44) for the symmetric case, i.e., a = 1, 
and the optimum asymmetry in the absence of fading as given by (41). 

Substituting a = 1 in (45), the parameters and. simplify to 


r = 4K 
4 1 + K + 4y 


= 


1 + K 

1 + K + 4y 


r = r = — 

s 2 ^3 1 + K + 2y 


1 + K 

1 + K + 2y 


(46) 


The curve labelled "symmetric" on Figure 7 is a plot of the upper bound of 
(44) combined with (46) as a function of the average bit energy-to-noise ratio 
E^/Nq with a Rician parameter K = 10 (typical of the mobile satellite channel). 

When the value of a in (41) is substituted in (45), after some 
simplification (44) and (45) can be written as 


P b± 


^1^3 exp [-(Cj + C 2 )] 
(1 - £2 ex P (-C 2 )) 


and 


I3;K 

r ' — . r ' 

M “ 1 + K + sj ’ ^i 


1 + K 


1 + K + ’ 1 


= 1, 2, 3 


. !b , S 

Q 1 = N Q ? °2 = an3 ? B 3 = N Q " an3 


(47) 


(48) 
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Figure 7. Bit Error Probability Performance vs E b /Ng f° r Rate l/ 2 
Trellis Coded QPSK in the Presence of Rician Fading; 

2 States, K = 10; Ideal Channel State Information 


The behavior of (47) combined with (48) is also illustrated Figure 7 by the 
curve labelled "optimum asymmetry (no fading)." 


For the Rayleigh case, the above results simplify even further. In 
particular, for the symmetric signal QPSK constellation, we get 


P b i 


1 + V 2N o 
V 2N o 


1 + W V + V 2N o, 


(49) 


whereas for the optimum asymmetry in the absence of fading. 


P b± 


1 + &n 3 
Q,n3 


1 + VV V 1 + V N o - an3 * 


(50) 
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These results are illustrated in Figure 8. Note that here the curve labelled 
"optimum asymmetry (no fading)" gives worse performance than that of the 
symmetric case. Thus, at least here, we clearly see the need for performing 
the asymmetry optimization in the presence of the fading. 

To determine the optimum value of a for the Rician case, we need to 
differentiate (44) with respect to a and equate the result to zero. This 
leads to a transcendental equation which must be solved numerically. Rather 
than doing that, it is more expedient to directly minimize (44) with respect 
to a using numerical techniques. When this is done, we obtain the optimum 
bit error probability bound labelled "optimum asymmetry (fading)" in 
Figure 7. For K = 10, this curve lies quite close to the "optimum asymmetry 
(no fading)" curve. However, as we can already deduce from Figure 8, this 
statement is not true for small values of K, in particular the Rayleigh 



Figure 8. Bit Error Probability Performance vs E^/Nq for Rate 1/2 
Trellis Coded QPSK in the Presence of Rayleigh Fading; 

2 States; Ideal Channel State Information 
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channel with K = 0. To exhibit the sensitivity of the optimum asymmetry 
condition to K, Figure 9 illustrates the optimum value of a as a function of 
E b /N 0 with K as a parameter. 

For the Rayleigh case, we can indeed determine the optimum asymmetry 
condition in closed form. In particular, differentiating (44) with K = 0 in 
(45) and equating the result to zero has the solution 


a = 


-* + Vo 


( v/l7 + 8^/Ng - L) 


4(1 ♦ Vo’ 


which when substituted back in (44) gives 


(51) 


(3 


+ J\7 + 8E l /NJ 


P b^ 


V"0' [4(E b /N 0> + (E b /N 0 


) 2 (7 + jl7 + 8E,_ /N_ ) J ‘ 


vo' 


16( v 17 + 8E U /N„ - 1)(E V /N n ) 2 (l + E b /N Q ) 4 


(52) 


b' 0 


b 0' 


This result is illustrated by the curve labelled "optimum asymmetry (fading)" 
in Figure 8 and is clearly superior to that corresponding to the symmetric 
signal constellation. 

Finally, for purpose of comparison, the corresponding upper bound on the 
performance of uncoded BPSK (same bandwidth as rate 1/2 trellis coded QPSK) in 
the presence of Rician and Rayleigh fading is also illustrated in Figures 7 
and 8. The analytical results for these curves are well known and are given by 

VVo\ 

l + K . E b /N 0 / 


P b i 


1 + K 

1 + K + E b / N() 


exp 
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V N 0' dB 

Figure 9. Optimum Asymmetry Parameter vs E^/Nq for Rate 1/2 

Trellis Coded QPSK in the Presence of Rician Fading 

for the Rician channel and 

p h < — : — < 54 > 

1 * V N 0 

for the Rayleigh channel. 

Depending on the shape of the bit error probability vs E^/Nq curve * 
one can often deduce some important practical conclusions by examining the 
asymptotic behavior of the curve. Since for the Rayleigh case (Figure 8), the 
error probability performance curves are essentially linear over a wide range 
of practical SNRs, one can approximately apply the asymptotic (large E^/Nq) 
result over this domain. In particular, for large E^/Ng, (49), (50), (52), 
and (54), respectively, become 
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< 


1 


(symmetric) 


(49’) 


'b‘ 


<^VV 


P ” " r ( _«tnj_ ) y 2 
L 1 + in 3 N oJ 


(optimum asymmetry - no fading) 


p < I 

” <W 


(optimum asymmetry - fading) 


and 


p. * 1 


b V N 0 


(uncoded) 


(50’) 


(52’) 


(54*) 


Thus, for example, comparing (49’) with (54’), we see that the effect of 
coding is to change the rate of descent of the error probability vs. E^/N^ 
performance from an inverse linear to an inverse square law behavior. If the 
QPSK constellation is now designed according to the optimum asymmetry for no 
fading, the performance is worse than that of the symmetric constellation by a 
factor of (1 + 8-n 3)/(/2 Q-n 3) or 1.3 dB. On the other hand, if the 
constellation is designed with the optimum asymmetry determined in the 
presence of fading, then, relative to the symmetric design, the performance 
is improved by a factor of Jz or 1.5 dB. From Figure 8, we see that these 
asymptotic results are almost achieved at an error rate of 10 


B. NO CHANNEL STATE INFORMATION 


Recalling (29b), we see that the transfer function T(D,I) of the 
pair-state transition diagram for the case of no channel state information is 

obtained by. replacing D^, this time by = D ^ £ n b ranc h 

label gains of (37a) where again 13 = 4, 4/(1 + a), or 4<r/(l + a) as appro- 
priate. Unfortunately, for the Rician distribution, the factor cannot 
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be evaluated in closed form. It can, however, be expressed as a single 
integral with finite limits as follows: 


D^Xfip _ e _ _1 I n (Q) exp [ n 2 (0) ] erfc n(®) d0] 


\S(E K /2N ) 

n(0) = - '/ic cos 0 

/l + K 


(55) 


This integral is easily evaluated using Gauss-Chebyshev techniques, i.e.. 


S' 


N 

n(0) exp tn 2 (0) ] erfc n(9) d0 = ~ n(9 k ) ex P erfc n(9 k ) 


k=l 


(56) 


where 0^ = (2k - l)ir/2N. 


For the Rayleigh case, we 
since, for K = 0, n(0) becomes 


can obtain a closed form result for this factor 
independent of 0. Thus, 


D 4X0P = 1 





(57) 


Evaluating T(D,I) of (37a) using (55) and performing the differentiation 
required in (35) gives the upper bound 


. -16X 2 (l+2a)/ (1+a) 

mm min 13 D 

X>0 a / ...2,,. v\2 

(l - ? 2 d‘ 16X /<Ua) ) 


(58) 
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where 


5. A e 
1 = 


-K 


;.af. 


(0) exp (n i (0)) erfc r^(0) d0 


/ir 0 

X J 3 i (E b /2N 0 ) 


n.(0) A 

Jl + K 


- 7k cos 0 


( 59 ) 


and B^; i = 1, 2, 3 are defined in (45). For the Rayleigh case one merely 
replaces ^ and rj^(0) of (59) by 


2 

A 1 - /rr n. exp (n.) erfc n- 
i = i r i i 

n i A XB i (E b /2N 0 ) 

and performs the same minimizations required in (58). 


(60) 


Figures 10 and 11 illustrate the analogous results to Figures 7 and 8 for 
the case where no channel state information is available.* Clearly, the lack 
of channel state information produces a noticeable degradation in system 
performance. To quantitatively assess this additional degradation (at least 
for the Rayleigh channel), we now derive asymptotic results analogous to 
(49'), (50'), and (52') for the no channel state information case. In 
particular, we use the asymptotic (large argument) expansion for erfc x, i.e.. 


erfc x ; sm <-y . ) 

/ir x 



in which case (57) simplifies to 


(61) 


Q 4XBp ~ 1 

2[XB(E b /2N Q )] 2 


(62) 


*For simplicity of presentation, we have chosen not to illustrate the results 
for the values of optimum asymmetry determined from the no fading analysis 
since we have already made the point that asymmetry should be optimized in 
the fading environment. 
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Figure 10. Bit Error Probability Performance vs Eb/Ng for Rate 1/2 
Trellis Coded QPSK in the Presence of Rician Fading; 

2 States, K = 10; No Channel State Information 


Using the appropriate values of B in (62), the branch gains of Figure 6 become 


a 


I 

4(2XE b /N 0 ) 2 


exp (4X 2 E b /Np) 



(63) 
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Figure 11. Bit Error Probability Performance vs E^/Ng for Rate 1/2- 
Trellis Coded QPSK in the Presence of Rayleigh Fading; 

2 States; No Channel State Information 


Substituting (63) into the transfer function of (37a) and performing the 
differentiation required in (55) gives the approximate upper bound on P 

b 

(valid for large E^/N^) 


P, < min 
b ~ X 


exp 


x V N o) 


4 <WV 4 < rh ) 


x 




exp 



(64) 
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Performing the minimization over X required in (64) (actually we minimize 
only the numerator since the denominator has little effect on this operation) 
gives 


opt 


2 (fvr)<VV 


which when substituted in (64) yields 


p b i 


{Hr 23 ) 


2 2 
e 


16 <VV 


Finally, the desired asymptotic results are 


V * 


9e 


16 W 


(symmetric) 


P, < 


2 V N o - *" 3 ' 


Vo - ln3 / 16 <E b /N o )2 


(optimum asymmetry - no fading) 




4 ( V N o v 


(optimum asymmetry - fading) 


( 65 ) 


( 66 ) 


(67) 
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SECTION VII 
SIMULATION RESULTS 


In this section, we describe and present the results of a software 
simulation of the system block diagram of Figure 2. For the application at 
hand, the development of a simulation has a manyfold purpose. First, it can 
be used to "verify" the theoretical results obtained in the previous section 
keeping in mind that the simulation is indicative of the exact system 
performance whereas the theoretical bit error rate expressions are upper 
bounds. Second, when the number of states in the trellis diagram becomes 
large (e.g., 16), determining the state transition diagram and its associated 
transfer function is a tedious task; in such cases, simulation is the more 
expedient approach. Finally, system degradation due to the finite size of 
interleaving and decoder buffer imposed by the practical constraint on the 
allowable total delay is analytically intractable, particularly when coupled 
with that caused by the "noisy" carrier demodulation reference produced by the 
pilot tone extractor. Hence, to predict true system performance corresponding 
to the real world environment, one must again turn to simulation. In the next 
paragraph, we expand upon the last of these issues. 

The block interleaver of Figure 2 can be regarded as a buffer with d rows 
which represent the depth of interleaving and s columns which represent the 
span of interleaving . Thus, the size of the interleaver (in symbols) is 
d x s. Data is written into the buffer in successive rows and read out of the 
buffer (the order in which it is transmitted over the channel) in columns. At 
the receiver, the block deinterleaver performs the reverse operation, i.e., 
the received soft quantized symbols are written into the buffer in successive 
columns and read out in rows. In practice, the interleaving depth should be 
chosen on the order of the maximum fade depth anticipated which, for the 
fading mobile satellite channel under investigation, depends on the doppler 
frequency or, equivalently, the vehicle speed. The smaller the doppler 
frequency, the longer the fade duration and vice versa. The interleaving span 
should be chosen on the order of the decoder buffer size. When this is done, 
the performance degradation (relative to that for the analytically tractable 
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assumption of infinite interleaving depth and buffer size) will be inversely 
proportional to the product of interleaving size and doppler frequency. 

On the other hand, when pilot tone(s) are used for coherent demodulation 
(as suggested in Figure 2), performance will degrade directly proportional to 
doppler frequency. The reason for this is that the bandpass filter(s) used in 
the pilot tone extractor to isolate the pilot tone(s) from the modulation must 
have bandwidth sufficiently wide to include the doppler shift. Thus, the 
larger the doppler, the wider the bandwidth of the filter(s) and hence the 
"noisier" the extracted demodulation reference. Assuming infinite 
interleaving and decoder buffer size, one can use the same analytical approach 
as previously discussed to derive upper bounds on the bit error probability in 
the presence of the noisy carrier reference. In particular, we first find the 
probability density function (p.d.f.) of the phase error associated with the 
demodulation reference signal which itself depends on the envelope of the 
fading sample. Next, the bit energy-to-noise ratio in the branch label gains 
(see Figure 6 for example) should be degraded by the cosine of the phase 
error, averaged over the phase error p.d.f., and then, as before, over the 
fading envelope p.d.f. Finally, using the averaged branch label gains in the 
transfer function bound T(D,I) gives the desired upper bound on average bit 
error rate. Even for the assumption of infinite interleaving and decoder 
buffer size, this computation is quite tedious. Thus, a simulation is 
preferable. 

Example 1 ; This example is for the "verification" of the analytical 
results obtained in the previous section. In particular, the system block 
diagram of Figure 2 was simulated for a rate 1/2, 2-state trellis code, 
symmetric or optimum asymmetric QPSK modulation, and either a Rician or 
Rayleigh fading channel. The interleaving size was chosen equal to 512 QPSK 
symbols (or equivalently, 512 input bits) which, for all practical purposes, 
approximates infinite interleaving. A doppler frequency of 100 Hz was chosen 
which makes the channel rapidly varying enough for the fading to be assumed 
independent from symbol to symbol. These two assumptions provide a memoryless 
channel as assumed in the analysis. Furthermore, the buffer size was chosen 
equal to 32 bits which approximates the assumption of an infinite bit buffer. 
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Numerical simulation results for this example are superimposed as dashed 
lines on the analytical results of Figures 7 and 8 corresponding respectively 
to the cases of Rician and Rayleigh fading with ideal channel state 
information. The discrepancy between the solid and dashed curves reflects the 
looseness of the upper bounds but the relative behavior of the analytical 
curves compares well with that of the simulation results. 

Example 2 : Here we consider the more practical case of a rate 2/3, 

16-state trellis code combined with symmetric BPSK modulation (it was shown in 
[5] that, for this case, the additional coding gain produced by the addition 
of asymmetry to the modulation is small and thus we have chosen to ignore 
it.) Although, with much computational effort and the assumption of infinite 
interleaving and buffer size, this system can be analyzed by the approach 
taken in the previous section, our interest here lies in computing the 
performance with limited interleaving and decoder buffer size as follows. 

At the present time, this system is a candidate for NASA's Mobile 
Satellite Experiment (MSAT-X) project [12] whose objective it is to transmit 
4800-9600 bps of digitally encoded speech over a 5-kHz RF channel with a bit 

_3 

error rate of 10 . To satisfy the previously mentioned constraint on total 

allowable delay, the interleaving size, interleaving depth, and decoder buffer 
size have been optimized at this bit error rate to achieve the minimum bit 
signal-to-noise ratio. For the specified delay constraint (60 ms), the size 
of the block interleaver and deinterleaver have been chosen equal to 128 BPSK 
symbols (or 256 input bits). With the above chosen interleaving size, the 
interleaving depth has been optimized by computer simulation and found to be 
equal to 16 symbols. Thus, the interleaving span is 128/16 = 8 symbols over 
the range of doppler frequencies from 20 Hz to 100 Hz (vehicle speeds of 
15 mph to 75 mph at UHF ) . Note, however, that for MSAT-X channels operating 
at low doppler frequencies such as 20 Hz, we can indeed have fade durations 
much longer than 16 symbols. In this, case, an interleaving size of 128 
symbols is not sufficient and, as we shall see shortly, a significant 
performance penalty occurs. Finally, with the above delay constraint imposed, 
the buffer size was optimized through simulation and found to be 32 symbols 
(or 64 bits). 
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Figures 12 and 13 illustrate the results of the simulation for perfect 

carrier and time synchronization and no intersymbol interferences (ISI). In 

particular. Figure 12 assumes a fixed doppler frequency of 100 Hz,' a fixed 

Rician parameter K = 10, and various interleaving and channel state 

information options. Also shown as a reference point is the performance of 

uncoded QPSK (same bandwidth as rate 2/3 trellis coded 8-PSK) at the chosen 
-3 

error rate of 10 . Figure 13 shows the effect of doppler frequency on 

system performance for the case of ideal channel state information and the 
same interleaving parameters as in Figure 12. 


Table 1 summarizes the above results by tabulating the required E^/N^ at 
a bit error rate of 10 ^ for each of the cases and also for Rician fading 
parameters of 5 and 7 dfi. From these numerical results, one can assess the 


< 

O' 

Q£ 

S 

Q£ 


| SIMULATION POINTS 


• RICIAN K = 10 dB 

• DOPPLER FREQUENCY = 100 Hz 

• INTERLEAVE = 128 SYMBOLS 

• DECODER BUFFER = 32 SYMBOLS 

• PERFECT SYNC (PHASE, TIME) 

• NO ISI 

• CSI = CHANNEL STATE INFORMATION 



UNCODED 

QPSK 


12 


BIT SNR 


Figure 12. Performance of 2/3, 16 State Code Over Fading Channel 
With/Without CSI With/Without Interleaving 
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BIT SNR 


Figure 13. Performance of 2/3, 16 State Trellis Coded 8-PSK Modulation 
(TCM) Over Rician Fading Channel With CSI and Interleaving 


coding gain (reduction in required E^/Ng relative to uncoded QPSK) 
achieved in each case. For example, without interleaving and channel state 
information, transmitting 2/3 coded 8-PSK over the K = 10 Rician channel 
produces a coding gain of 1.6 dB. When 128 symbol interleaving is added, the 
coding gain is increased to 3.1 dB. If, in addition to interleaving, we 
provide ideal channel state information, another 0.5 dB coding gain is 
achieved which brings the total coding gain to 3.6 dB in this particular 
environment. As is obvious from Figure 12, this coding gain would be greater 
at lower bit error rates. When the doppler is decreased from 100 Hz to 20 Hz, 
Figure 13 and Table 1 show a reduction in coding gain of 1 dB due to the 
limitation on the size of the interleaver. 
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Table 1. Summary of Results 


Type of 
Modulation 

Fading 
K, dB 

Doppler 
K, dB 

Block 
Interleave 
128 Symbols 

Channe 1 
State 

Information 

(CSI) 

Required Bit 
SNR, dB at 
BER = 10“ 3 

QPSK 

10 

100 

No 

No 

9.6 

TCM 

10 

100 

No 

No 

8.6 

TCM 

10 

100 

Yes 

No 

6.6 

TCM 

10 

100 

Yes 

Yes 

6.0 

TCM 

10 

20 

Yes 

Yes 

7.0 

QPSK 

7 

100 

No 

No 

12.0 

TCM 

7 

100 

Yes 

Yes 

7.0 

TCM 

7 

20 

Yes 

Yes 

8.5 

QPSK 

5 

100 

No 

No 

15.0 

TCM 

5 

100 

Yes 

Yes 

7.8 

TCM 

5 

20 

Yes 

Yes 

10.5 

*Decoder 

buffer size 

= 32 symbols. 




As previously mentioned, all of the results in Figures 12 and 13 and 
Table 1 are for the case of perfect carrier synchronization. Using a dual 
tone calibration technique (DTCT) [13] wherein two tones of equal power are 
inserted symmetrically at the edges of the data spectrum for the purpose of 
coherent demodulation, we have found by simulation that the noisy carrier 
reference produced by the appropriate pilot tone extractor [13] produces about 
a 2- to 3-dB degradation in performance depending on the value of the Rician 
parameter K. This is caused by the fact that the bandwidth of the pilot tone 
bandpass filters has to be chosen wide enough to accommodate the maximum 
doppler of 100 Hz and the total power of the two pilots is 7 dB below that of 
the data (this ratio has been shown to be optimum). Simulation results with 
the DTCT technique are shown in Figure 14 for Rayleigh and various Rician 
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Figure 14. Performance of TCM vs QPSK With DTCT 



channels under the assumption of ideal channel state information, 128 symbol 
interleaving, and a doppler frequency of 20 Hz (worst case for our applica- 
tion). Note that in this simulation, the bandwidth of the pilot bandpass 
filters was chosen to accommodate the 100-Hz doppler as would be the case in 
the actual hardware design. 
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